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Abstract 

We consider a recent proposal by Horvath et al. to address the question 
whether topological charge fluctuations in QCD are instanton dominated via 
the response of fermions using lattice fermions with exact chiral symmetry, 
the overlap fermions. Considering several volumes and lattice spacings we 
find strong evidence for chirality of a finite density of low-lying eigenvectors 
of the overlap-Dirac operator in the regions where these modes are peaked. 
This result suggests instanton dominance of topological charge fluctuations in 
quenched QCD. 
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I. INTRODUCTION 



Understanding the mechanisms of confinement and chiral symmetry breaking in QCD 
has been a major goal of non-perturbative field theory for many years. A good candidate 
for the source of chiral symmetry breaking are topological fluctuations of the gauge fields 
associated with instantons and anti-instantons. Through mixing, the "would be zero modes" 
from the instantons and anti-instantons become the small modes that lead, when present 
with a finite density, to chiral symmetry breaking. Many of the low energy properties of QCD 
can be explained phenomenologically by the interactions of the fermions with instantons and 
anti-instantons |TJ,|2|]. To our knowledge, confinement, though, is not one of those properties. 

However, as pointed out by Witten || , there is an inconsistency between instanton based 
phenomenology and large-iV c QCD. Instantons would produce an rj' mass that vanishes 
exponentially for large N c , while considerations based on large- N c chiral dynamics suggest 
that the rf mass should be of order 1/N C . The topological charge fluctuations then should be 
associated not to instantons but rather some other, confinement-related vacuum fluctuations. 
And the strong attraction produced by these confinement-related fluctuations would also 
induce the breaking of chiral symmetry. 

In an attempt to settle the issue, whether topological charge fluctuations are instanton 
dominated or not, Horvath, Isgur, McCune and Thacker [Q investigated the response of 
fermions by considering low lying eigenmodes of a lattice discretized Dirac operator and 
quantifying the extend to which, in the peak regions of these eigenmodes they are chiral, 
as would be expected if the almost zero modes around instantons dominate these low lying 
eigenmodes. They introduced a local "chiral orientation" parameter, or "chirality" param- 
eter, X(x), via the definition 



where and iPr{x) are the left- and right-handed components of the eigenmode. Chiral- 

ity near ±1 would be obtained from almost zero modes due to instantons and anti-instantons, 
lifted from being exact zero modes by their mixing. 

Horvath et al. Jlj chose the Wilson discretization for their lattice Dirac operator and 
claim to have found evidence against instanton dominance of topological charge fluctuations 
from the behavior of the chirality parameter. However, the Wilson-Dirac operator explicitly 
breaks chiral symmetry at any finite lattice spacing a. As a consequence, the Wilson-Dirac 
operator does not have exact (chiral) zero modes in gauge fields with a non-trivial topological 
charge. And the discretization effects responsible might well also strongly affect the chirality 
parameter. To our knowledge, the magnitude of such distortion effects at the presently 
accessible lattice spacings are not known. However, the scaling violations for hadron masses 
due to the 0(a) lattice artifacts at the lattice spacing used by Horvath et al. (a ~ 0.17 fm 
at (3 — 5.7) are known to be as large as 30-40% ||. 

In a response to M DeGrand and Hasenfratz presented evidence that is suggestive of 
instanton dominance of topological charge fluctuations ||. They used a variant of the 
overlap-Dirac operator, proposed and investigated in [[7], for the lattice Dirac operator. The 
overlap-Dirac operator ||, and any variant that substitutes the Wilson-Dirac operator used 
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in its construction by another suitable Wilson-like discretization of the Dirac operator, has 
the chiral and topological properties of the continuum Dirac operator even at finite lattice 
spacing, e.g. exact chiral zero modes in gauge fields with a non-trivial topological charge. 
However, DeGrand and Hasenfratz coupled their overlap fermions not directly to the "rough" 
lattice gauge fields, but rather to "APE smeared" gauge fields, essentially Gaussian smeared 
fields with fixed width in lattice units. While such a fermion action is local and therefore 
does not change universality class, the distortions the smearing procedure might induce at 
finite lattice spacings are not known. Indeed, the level of smearing used by DeGrand and 
Hasenfratz removes ultraviolet fluctuations in the gauge field to an extent where pure gauge 
observables start to "feel" and be able to identify instanton-like topological excitations that 
on the "rough" lattice gauge fields are completely obscured by the dominant ultraviolet 
fluctuations. For example, the naive topological charge measurements gives results that are 
clearly peaked around integer values. The smearing used distorts the short distance part of 
the heavy quark potential, to a distance of about 3 lattice spacings, and lowers the value 
obtained for fg when the fat links are used in the clover Dirac operator ||. Whether the 
smearing affects properties of light fermions has, to our knowledge, not been investigated 
thoroughly. 

In an earlier paper, using the same variant of the overlap-Dirac operator, DeGrand and 
Hasenfratz produced evidence that the chiral density u(x) = ip^^ip^x) comes in lumps on 
the lattice , and furthermore that the correlation function of the chiral density with itself 
is similar in shape to the correlation function of the chiral density u{x) with the topological 
charge density Q(x) obtained from an APE-smeared operator. But, since APE-smearing 
is involved, these results are not beyond doubt. Indeed, the authors of Ref. || argue that 
the "lumpiness" observed in the correlation function of the chiral density uo{x) could occur 
without these lumps being purely right-handed or left-handed. 

Here we address these questions using the standard overlap-Dirac operator || with the 
fermions coupled to the original lattice gauge fields. We thereby avoid potential distortions 
due to smearing while we retain all chiral properties including existence of exact zero modes 
and thus avoid the lattice artifacts, unknown in magnitude and significance, of Wilson 
fermions. 



II. SETUP 

We consider the standard massless overlap-Dirac operator M 

D(0) = ±[l + l5 e(H w (M))] (2) 

with H W (M) = jD w (-M) and D W (M) the usual Wilson-Dirac operator. Then, if 2 (0) = 
D*(0)D(0) commutes with 75 and can be simultaneously diagonalized [11 . 1^] . H 2 (0) can 



have zero eigenvalues with chiral eigenmodes, which are also eigenmodes of -D(O), due to 
global topology. The non-zero eigenmodes of H 2 (0) are doubly degenerate and have opposite 
chirality, 

# 2 (°)^U = *Vt4 > 7s^T = V>t >7sVn = > ( 3 ) 
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with < A < 1. The non-zero (right) eigenmodes of -D(O) are then easily obtained as 

V>± = 4=(Vt±#|) (4) 



with eigenvalues 



A± = A 2 ± iXVl^X 2 . (5) 
The "chirality" parameter of eq. ([I]) is then given by 



tan(^(l + X(x))) = ^ mK ' , (6) 
V4 V n ) \^Jx)^(x)J 



and is equal for the two related modes ip±. In the following they will be therefore counted 
as one mode. The exact zero modes, of course, have X(x) = ±1. 

We computed low-lying eigenmodes ip^^ of H 2 (0) with the Ritz functional algorithm of 
Ref. |13| . For the sign function e(Hw) in eq. (H) we used the optimal rational approximation 



of [|TT|,|l3l with "projection of low-lying eigenvectors of H w " to ensure sufficient accuracy. 

III. RESULTS 

We computed the 20 eigenmodes if>^^ of H 2 (0) with smallest A on ensembles of 8 3 x 16 
pure gauge Wilson action configurations with f3 = 5.7 and 5.85, and Wilson-Dirac mass 
of M = 1.65. The gauge configurations have lattice spacing of about 0.17 and 0.125 fm, 
respectively, and hence volumes of about 7 and 2 fm 4 . We also computed the 20 lowest 
eigenmodes on an ensemble of 6 3 x 12 lattices with (3 = 5.7, with a volume of about 2.2 fm 4 , 
fairly close to the volume of the (3 = 5.85 ensemble. 

In Fig. |1] we show the chirality histograms for the lowest two non-zero modes at the 
2.5% of the lattice sites with the largest ip^ip(x). In this, and in all similar figures, the 
histograms are normalized such that the area under the histogram is equal to the fraction of 
sites considered. In their peaks, the lowest two eigenmodes are fairly chiral. Comparing the 
two ensembles with the same gauge coupling (and hence lattice spacing) we see a dramatic 
dependence on the physical volume. In an instanton liquid picture of the vacuum, the number 
of instantons and anti-instantons, and hence the number of almost zero modes grows linearly 
with the volume. Then it is not surprising that the lowest modes become increasingly chiral 
in their peaks with increasing volume. Comparing the two ensembles with almost equal 
volume in physical units, we notice a tendency of the peaks to become more chiral as the 
lattice spacing is decreased. 

Comparing with Fig. 1 of Ref. M we see that in their peaks the lowest two non-zero 
eigenmodes of the overlap-Dirac operator appear to be more chiral than the real modes of 
the Wilson-Dirac operator. Considering that Ref. || used a lattice with larger volume of 
about 35 fm 4 and the dramatic volume dependence of the chirality histogram observed here 
leads one to suspect that the chiral symmetry breaking lattice artifacts of the Wilson-Dirac 
operator are truly significant at the lattice spacing used. 
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FIG. 1. Chirality histograms for the lowest two non-zero modes of the overlap-Dirac operator 
at the 2.5% sites with the largest ip^ip{ x ) on the three ensembles with Wilson gauge action. 
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FIG. 2. Chirality histogram similar to Fig. |l|, but for the lowest four non-zero modes at the 
0.5% sites with the largest ip'ip(x). 

Comparing with Fig. 1 of Ref. ||, on the other hand, where a physical volume of about 
3.5 fm 4 was used, indicates that the APE smeared fields used by DeGrand and Hasenfratz 
enhance the chirality in the peaks of the lowest-lying non-zero modes somewhat. DeGrand 
and Hasenfratz already noted a slight dependence on the smearing of gauge fields (see Fig. 4 
of H), though the effect appears much smaller than the effects of the lattice artifacts from 
using Wilson fermions. 

Considering also the next two higher modes while further restricting the peak region 
gives the chirality histograms for the two smaller lattices shown in Fig. These modes 
appear still rather chiral in their peaks. On the larger volume, as expected, a larger number 
of modes are chiral in their peak, as seen in Fig. 

To confirm the trends with volume and lattice spacing, we would like to consider a still 
larger volume, as well as a lattice at smaller lattice spacing. Unfortunately, the numerical 
implementation of overlap fermions becomes increasingly difficult and costly, due to the finite 



density of modes with almost zero eigenvalues of the underlying Wilson-Dirac operator |L4 



It is believed that these modes are caused by lattice artifact dislocations in the gauge fields 



generated with Wilson like gauge actions. It is known [|i5H l7|l that these same small modes 
are responsible for the residual chiral symmetry breaking observed for domain wall fermions 
with finite 5-th dimension L s . It was found that improving the gauge action via the Symanzik 
improved action suppresses the dislocations and decreases the density of small modes |TJ| . It 
has also been observed that the residual chiral symmetry breaking for domain wall fermions 



on gauge fields generated with an improved gauge action is much reduced [|T^[IB]|n| . And 



in Ref. |20[ it was noted that the cost of implementing the overlap-Dirac operator on gauge 
fields with a tree-level Symanzik improved gauge action is lowered. 

We made a study of various gauge actions to find which one lowered the density of small 
modes and would consequently lower the cost of implementing the overlap-Dirac operator. 
The density of zero eigenvalues of the hermitian Wilson-Dirac operator p(0) at the mass M = 
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FIG. 3. Chirality histogram for the lowest six non-zero modes at the 2.5% sites with the largest 
i/j^ip(x) on the lattice at j3 = 5.7 with the larger volume. 



1.65 was determined for quenched backgrounds from the Iwasaki [22], tree-level tadpole- 



improved and 1-loop tadpole- improved Liischer-Weisz gauge action [23]. This was done 



by computing the low lying eigenvalues of H W (M) using the Ritz functional ]T3|]. The 



low lying eigenvalues over the whole ensemble are then used to obtain the integral of the 
spectral density function, J Q p(X')dX'. A linear fit in A is made, and p(0) is obtained as the 
coefficient of the linear term. Comparisons were made with the Wilson gauge action at two 
scales determined by the string tension. The results shown in Table | show the quenched 
Iwasaki gauge action to have the lowest p(0). For example, at the lattice spacing equivalent 
to Wilson gauge action (3 = 5.7 the p(0) is reduced by a factor of about 5. 

We therefore generated gauge field ensembles with Iwasaki's action with (3 = 2.2872 and 
2.45, chosen to give the same lattice spacing, in units of the string tension, as the ensembles 
already considered with Wilson action at (3 = 5.7 and 5.85. As expected, since the density 
of small modes of the Wilson-Dirac operator used in the overlap fermion construction were 
considerably reduced, the overlap fermion computation became much faster. For example, at 
the lattice spacing equivalent to Wilson /? = 5.7 and lattice size 8 3 x 16 the cost of the inner 
conjugate- gradient step in applying the overlap-Dirac operator was reduced from between 
200 to 400 iterations to between 120 to 170 iterations. The average time per configuration 
to compute the 20 lowest eigenvectors was reduced by roughly a factor of 3. This made it 
feasible to consider lattices of size 12 3 x 16, giving a volume of about 23.6 fm 4 at (3 = 2.2872 
and 6.7 fm 4 at (3 = 2.45, i.e. about the same volume as the 8 3 x 16 lattice at (3 = 2.2872. 
We also generated a 12 3 x 16 at (3 = 2.65, chosen such that the volume is about 2.1 fm 4 , 
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TABLE I. Comparison of the density of zero eigenvalues p(0) of the hermitian Wilson-Dirac 
operator, in lattice units, for various gauge actions. The comparisons were made at the same 
lattice spacings as the Wilson gauge action (3 = 5.85 and 5.7. The scale was set using the string 
tension. Shown are the results for the Wilson, Iwasaki, tree-level and 1-loop tadpole-improved 
Liischer-Weisz gauge actions. The string tensions for the Wilson gauge action were taken from 
Ref. ||21j] . The string tension for the 1-loop tadpole- improved Liischer-Weisz gauge action comes 
from an interpolation of previous results, while the others were computed in this work. All tests 
used a Wilson-Dirac mass of M = 1.65. The Iwasaki gauge action has the lowest p(0) at fixed 
a^/a. 



action 


P 


V a 4 




p(0) x 10 5 


Wilson 


5.70 


8 3 x 16 


0.3917 


162(9) 




5.85 


8 3 x 16 


0.2864 


21(3) 




6.00 


24 3 x 32 


0.2197 


2.3(2) 


Iwasaki 


2.2872 


8 3 x 16 


0.3885(93) 


32(4) 




2.45 


12 3 x 16 


0.2840(53) 


2.4(9) 


0-loop TI LW 


7.26 


8 3 x 16 


0.3900(76) 


58(6) 




7.60 


12 3 x 16 


0.2872(39) 


4.3(8) 


1-loop TI LW 


7.79 


8 3 x 16 


0.3917 


44(5) 




8.13 


12 3 x 16 


0.2864 


3.9(7) 



the same volume as the smallest lattice at the other two gauge couplings. The lattice size 
and volume, both in lattice and physical units, as well as the average number of zero modes 
and the fluctuation of the topological charge, obtained from the number of zero modes, of 
all ensembles considered are listed in Table O- 

The chirality histograms of the two lowest non-zero modes at the 2.5% of sites with largest 
ip^ip{x) for all six ensemble with Iwasaki glue are shown in Fig. |4j. The physical volume of 
the systems shown in the left column are approximately the same. Comparing the three 
histograms confirms the trend, already observed with Wilson glue, that the lowest non- 
zero modes tend to become more chiral in their peaks as the lattice spacing is decreased. 
This strongly suggests that the behavior observed here will survive the continuum limit. 
Comparing the three histograms for j3 = 2.2872, on the other hand, confirms the strong 
volume dependence already observed with Wilson glue. 

In Fig. |S| we show the chirality histogram of the lowest six non-zero modes on the two 
systems with physical volume of about 7 fm 4 and in Fig. |6| the chirality histogram from all 
non-zero modes out of the 20 lowest modes that we computed (there are between 11 and 20 
non-zero modes per configuration). Comparing with the histograms in the left column of 
Fig. f| indicates that the number of non-zero modes with similar chirality histograms grows 
roughly like the physical volume so that there is a finite density of modes which are chiral 
in their peaks. This is good evidence that we are not observing a finite volume effect that 
would disappear in the infinite volume limit. 

In the chirality histograms shown we typically considered the contribution from the 2.5% 
sites with the largest ijj^ip{x). In Table |T1| we list the sum of ip^ip{x) over those sites as a 
percentage of the sum over all sites for some low-lying eigenvectors for all the gauge field 
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TABLE II. The gauge ensembles with Wilson glue in the first three lines and with Iwasaki glue 
in the rest. Listed are the volume in lattice units and in units of fm 4 , the number of configurations 
analyzed, the average number of zero modes and the average of the square of the global topological 
charge as obtained from the number of zero modes. 






V [a 4 ] 


V [fm 4 ] 


N cf 9 


<IQI> 


(Q 2 ) 


5.70 


6 3 x 12 


2.2 


50 


1.2(2) 


2.8(5) 


5.70 


8 3 x 16 


7.0 


48 


2.0(2) 


7(1) 


5.85 


8 3 x 16 


2.0 


93 


1.3(1) 


2.6(4) 


2.2872 


6 3 x 12 


2.2 


100 


1.0(1) 


1.7(2) 


2.2872 


8 3 x 16 


7.0 


69 


2.0(2) 


5(1) 


2.2872 


12 3 x 16 


23.6 


10 


4.2(9) 


25(8) 


2.45 


8 3 x 16 


2.0 


197 


1.0(1) 


1.8(2) 


2.45 


12 3 x 16 


6.7 


47 


2.3(4) 


8(1) 


2.65 


12 3 x 16 


2.1 


50 


1.3(2) 


2.4(6) 



TABLE III. The average of the sum of ip'ip(x) over the 2.5% of sites kept for the chirality 
histograms as a percentage of the sum over all sites for some low-lying eigenvectors for the various 
gauge field ensembles. 






V[a*] 


ev 1 


ev 2 


ev 4 


ev 6 


ev 20 


5.70 


6 3 x 12 


8.7(3) 


7.9(3) 


7.0(2) 


6.4(2) 


5.1(1) 


5.70 


8 3 x 16 


11.4(3) 


11.0(3) 


9.7(2) 


9.3(2) 


6.9(5) 


5.85 


8 3 x 16 


10.6(3) 


9.5(3) 


7.9(2) 


7.5(2) 


6.1(1) 


2.2872 


6 3 x 12 


8.9(2) 


7.8(2) 


6.9(2) 


6.3(1) 


5.1(1) 


2.2872 


8 3 x 16 


11.5(3) 


10.6(2) 


9.3(2) 


8.5(2) 


6.2(2) 


2.2872 


12 3 x 16 


12.9(5) 


13.4(6) 


12.3(4) 


12.0(4) 


11.3(5) 


2.45 


8 3 x 16 


9.9(2) 


8.7(1) 


7.7(1) 


7.2(1) 


6.0(1) 


2.45 


12 3 x 16 


12.4(3) 


11.5(3) 


10.3(2) 


9.9(3) 


7.8(3) 


2.65 


12 3 x 16 


10.4(3) 


9.2(2) 


8.0(1) 


7.7(1) 


6.9(3) 
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FIG. 4. Chirality histograms for the lowest two non-zero modes of the overlap-Dirac operator 
at the 2.5% sites with the largest ip^ip(x) on the six ensembles with Iwasaki gauge action. The 
systems in the left column all have approximately the same volume in physical units. The systems 
in the top two panels in the right column also have the same, about a factor of 3 larger, volume, 
while the system in the lower right hand corner panel has the largest volume. 



10 



0.03 



0.02 



0.01 - 



0.00 



(8=2.2872, 8 3 xl6 




0.03 



0.02 



0.01 - 



0.00 



1 1 1 I 1 1 
(9 = 2.45, 12 3 xl6 






X 



FIG. 5. Chirality histogram similar for the lowest six non-zero modes at the 2.5% sites with 
the largest ip^ip(x) on the two lattices with volume of about 7 fm 4 . 
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FIG. 6. Chirality histogram for all the non-zero modes out of the 20 lowest modes at the 2.5% 
sites with the largest ^ij){x) on the lattice at (3 = 2.2872 with the largest volume. 
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FIG. 7. Chirality histogram for the lowest ten non-zero modes of the overlap-Dirac operator at 
the 6% sites with the largest ip^ij)(x) on 24 2 U(l) configurations at j3 = 1.9894. 

ensembles considered. The more chiral the eigenvectors are on those 2.5% sites, the higher 
is the percentage of the eigenvector contained on those sites. 

IV. ABELIAN GAUGE THEORIES 



To contrast the results for quenched QCD of the previous section we considered the 
chirality parameter also for quenched U(l) theories in two and four dimensions. In 2d, also 
considered in Ref. Q], topology plays a crucial role and one expects local chirality in the 
peaks of low-lying modes. This is clearly seen in Fig. [7], where the 10 lowest non-zero modes 
have been kept. Comparing with Fig.s 6 and 7 of Ref. [|J demonstrates again the superiority 
of the overlap-Dirac operator over the Wilson-Dirac operator for the purpose of investigating 
chirality properties. The histogram in Fig. [7| for the 10 lowest non-zero modes shows better 
chirality peaks than Fig. 6 of [|J for the real modes of the Wilson-Dirac operator, the "should 
be zero modes" , were it not for the explicit chiral symmetry breaking from the Wilson term. 

In 4d U(l), in the confined phase, chiral symmetry is spontaneously broken, and a finite 
density of near zero modes exists. But there are no instantons, and so one would not 
expect the near zero modes to be dramatically chiral even in their peak regions. There 
do exist exact zero modes of the overlap-Dirac operator, which again, of course, are chiral, 
although their origin is not completely understood |23]]. We analyzed, stored eigenmodes 
from Ref. | 24]| . The chirality histogram for the lowest two non-zero modes is shown in Fig. || 



We note that while there is some mild indication of chirality peaking, the proportion of sites 
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FIG. 8. Chirality histogram for the lowest two non-zero modes of the overlap-Dirac operator at 
the 1% (left panel) and 0.25% (right panel) sites with the largest ip^ip(x) on 8 4 U(l) configurations 
at = 0.9. 

showing this behavior is much reduced compared to the SU(3) case and appears to be of a 
qualitatively different nature than seen before. Also, the sum of ip^ip{x) over the sites kept 
for the histograms of Fig. || where almost identical for all 12 lowest eigenmodes (4.0% and 
1.20%, respectively; for the 2.5% of sites with largest ip^ip(x) the sum is about 8.4% for all 
eigenmodes as compared to the non-abelian case in Table [TTJ). m the U(l) case the (near) 
lack of peaking could conceivably come about from the scenario outlined by Horvath, et. al. 
- namely from the confinement inducing vacuum fluctuations. 



V. CONCLUSIONS 

Horvath et al. proposed an observable, the "chirality" parameter, designed to answer 
the question whether topological charge fluctuations in (quenched) SU(3) gauge theory are 
instanton dominated or not by the response of fermions. Chirality of the low lying modes near 
±1 in the regions where they have the largest magnitude would signal instanton dominance. 
Using Wilson fermions, afflicted by potentially strong lattice artifacts, they claimed evidence 
that the topological charge fluctuations are not instanton dominated. We took up this 
question using overlap fermions, a lattice fermion formulation known to have the same 
chiral properties as continuum fermions, and therefore much less afflicted by lattice artifacts 
and better suited to the study of properties closely related to chiral symmetry. 

Considering several volumes and lattice spacings, we found convincing evidence for chi- 
rality of the low lying modes in their peak region, confirming results of Ref. || obtained 
using overlap fermions coupled to smoothened gauge fields. We note that improving the 
gauge action did not significantly change the degree of chirality at fixed lattice spacing and 
lattice size. Our results give evidence that the number of modes which are chiral in their 
peaks grows linearly with the volume so that there is a finite density of such modes, and 
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that the chirality of the modes becomes more pronounced as the lattice spacing is decreased. 
Therefore, our observations should remain valid in the continuum limit. 

That the results of Ref. M are strongly affected by lattice artifacts was confirmed by 
the very recent paper of Hip et al. |[25|| . These authors consider an improved chirality 
parameter for Wilson fermions, clover improved Wilson fermions, and lattices at smaller 
lattice spacing than Ref. Q . With these reductions of lattice artifacts, they conclude that 
instanton dominance of topological charge fluctuations is not ruled out by the response of 
(improved) Wilson fermions. Together with the results from overlap fermions presented in 
this paper and the lack of significant chirality enhancement in the 4-d U(l) model where 
instantons should not exist, the case for instanton domination in 4-d SU(3) gauge theory, 
as measured by the chirality parameter, becomes even more compelling. 
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